Non-mean-field effects in systems with long-range forces in competition 
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We investigate the canonical equilibrium of systems with long-range forces in competition. These 
forces create a modulation in the interaction potential and modulated phases appear at the system 
scale. The structure of these phases differentiate this system from monotonic potentials, where only 
the mean-field and disordered phases exist. With increasing temperature, the system switches from 
one ordered phase to another through a first-order phase transition. Both mean-field and modulated 
phases may be stable, even at zero temperature, and the long-range nature of the interaction will 
lead to metastability characterized by extremely long time scales. 



Long-range interactions are present in a wide variety 
of physical systems such as gravitation, plasma physics, 
wave-particle interactions or two-dimensional fluids 
The dynamics and thermodynamics of these systems is 
dominated by a mean-field phase where all elements in 
the system are in a similar phase (e.g. polarization of 
spins). In lattices with long-range forces,this result was 
confirmed in microcanonical simulations ||, and rigor- 
ously for microcanonical and canonical ensembles [JQ. 
These results validate the mean-field approach for long- 
range systems. 

The scenario changes when a short-range force enters 
into competition with the long-range one. It was shown 
that stripped patterns and other modulated phases may 
emerge |6l— ll 31] . A similar scenario occurs in the case 
of competing short-range forces: superconducting films, 
magnetic systems, chemical reactions (i.e. Turing pat- 
terns) and copolymers 1J] . Short-range purely repulsive 
forces can also generate such a phenomenon |l5| . 

Systems with two competing long-range forces are in- 
vestigated here, producing a situation where phases mod- 
ulated at the size of the system appear. This study is 
motivated by proposals of tablctop experiments where 
long-range forces are created artifically using cold atoms 
or colloids (l6l - [l8j . For example, Chalony et al. recently 
took advantage of laser- induced atom-atom interaction to 
create an artificial one-dimensional gravity These 
works open new perspectives to study, and tune, long- 
range forces. Modulated long-range potentials are also 
naturally found in the interaction betweem light and ul- 
tracold matter [2(J 21 1, making these systems appealing 



for future work. 

In this paper, it is shown that the competition of long- 
range forces in a lattice leads to canonical equilibria with 
phases modulated at the scale of the system. The free 
energy landscape is strongly modified and exhibits sev- 
eral minima for different ordered phases. As a result, 
the equilibrium phase depends on the temperature and 



phase transitions between macrophascs appear. Addi- 
tionally, minima that do not correspond to the canonical 
equilibrium give rise to metastable states characterized 
by extremely long lifetimes. 

To investigate this force competition, an Ising-like 
model of A spins on a lattice is used, inspired from the 
Dyson model with a single force |22j-|24|. It is a simple 
long-range model that captures the main features of force 
competition, and it is represented by the Hamiltonian 



(1) 



Each spin can take values <x; = ±1 and is at a position 
ri = i/N on a lattice. The coupling between two spins, 
Jij is a function of their mutual distance, and can be 
characterized by a competition of power laws 



J l} = V(dij) = 



\i-3\\ ai 



\i - j\\ a 2 



(2) 



Periodic boundary conditions are considered ||i — j\\ = 
, 1 — |rj — Tj |), though the general conclusions 
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of this work should remain valid for other boundary con- 
ditions. Using the eigenbasis on the potential (i.e. the 
Fourier basis) 13| the Hamiltonian ([I]) reads 



N 



(3) 



where the eigenmodes and eigenvalues are 

o c ,k = N^ 1 J2iCos{2Trkr i )a i , 
CT s ,fc = A -1 ^iSin(27rfcri)(7i, 
Afe = A -1 J2i cos (27rfcr-i)Ji i j 



(4) 



A single power-law potential yields a monotically de- 
creasing spectrum, but two competing ones can break 
this monotonicity, and allow for the emergence of mod- 
ulated phases [13[ shown in FigJT] The partition func- 
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FIG. 1. Spectrum of a single power-law potential (black plus), 
with a = 0.7, and of the two competing power-law potentials 
under study (blue crosses), with a = 0.1 and a' — 0. 
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tion, Z = ,_ ±1 exp(— /3if), contains a quadradic de- 
pendence on a in the exponential, which can be turned 
into a linear one using the Hubbard-Stratanovitch trick 



exp 



fNf3\d< 



N(3\ 



2tt 



dy exp 



N(3Xy 2 



Nl3Xy3j 
(5) 

where /3 ~ T 1 . Considering the K first Fourier modes, 
with K a cut-off on the Fourier modes at an arbitrary 
order, and applying Eq.Q to them 



1 (Nl 



if/2 



^ / dy k dy' k exp 

LI J 



<Ti=±l 



. (6) 



with A 



Afe. Polar coordinates, (yk = 
Pk cosyjfc, y' k = pkshnpk), are introduced and the sum 
over the configurations of spins is performed 

I (N£\ K ' 2 



Z = 



A V 2tt J 



log cosh 



PkdpkAkpk exp 



■(7) 



ft is noteworthy that the quantity in square bracket 
is extensive, so Z can be written as Z oc 
J dV exp[— N(f>(/3, pk-, <pk)] where 4> & function that is well- 
defined in the thermodynamic limit 



fll/2 l0 § C0Sn (Z 3 Efe A fe/°/c cos ( 27rfcr 



Vk))dr. 
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A similar treatment can be performed for short-range 
interactions using a mean-field approximation, but in the 
long-range case, it is an exact result. Using to the saddle- 
point method, in the thermodynamic limit, the canonical 
equilibrium is given by the minima of the function <j) 



Pk 



f 1/2 ( 

/ dr cos (2irkr) tanh (3 \, ^fe' Pk' cos (27r/c'r) 

J-X/2 V k' 



for each mode k. We have here omitted the dependence 
on ifk, since the equilibrium is degenerate on this vari- 
able. Indeed, for periodic boundary conditions, the par- 
ticles are not distinguishable and the equilibrium has a 
translational degeneracy. 

In general, Eq.© yields solutions that arc a combi- 
nation of eigcnmodes. For the sake of simplicity, let us 
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first consider a potential with only two modes k = 0, 1. 
With the eigenvalues restricted by Afc>2 = 0, the min- 
ima are determined by Eq.®. For k = 0, 1, the set of 
two equations is closed, therefore the amplitude of the 
k > 2 modes may be derived from these. The function 
has two minima, p and ~p 1 , which are associated to each 
mode 

p Q = tanh(/3A p ) , (9a) 

,1/2 

Pi = I dr cos (2-kt) tanh (fi\ip l cos (27rr)) . (9b) 

7-1/2 

In particular, using Eq.Q, one can show that the equi- 
librium value of eigenmode k is given by: 



2 dlogZ __ 2 



(10) 



Here, ~p k must correspond to the minimum of the canon- 
ical equilibrium. 

In the case of a single power-law potential, V(r) ~ 
l/r a , the series of the A& monotically decreases and 
the absolute minimum is associated to the mean-field 
mode, differing from above the critical temperature 
where the system remains in a disordered phase (i.e. 
Pk = 0, Vfc) Q. However, when power-law forces are 
competing, provided Ai > 2Ao, the absolute minimum 
yielded by © depends on the temperature. In fact, 
the critical temperature below which a mode can exist 
is given by T k = Afc(l + 5o,fc)/2. Hence, if Ai > 2A , the 
first Fourier mode exists at temperatures T§ < T < Tf 
where the mean-field mode does not. 

This phenomenon is illustrated in Fig 121 where the pro- 
file of 4> is plotted. Above the critical temperature Tf, 
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only the minimum corresponding to the disordered phase 
Po = Pi = exists, Figj2ja). Below Tf , but above T§, a 
minimum associated to the first Fourier mode ~p 1 > ap- 
pears and represents the statistical equilibrium, Fig(2Jb) . 
When the temperature reaches Tq , a minimum associated 
to the mean-field phase, ~p > and ~p 1 = 0, appears, but 
is still higher than the one associated with k = 1 Fourier 
mode, while the latter remains the canonical equilibrium, 
Fig Etc). However, at some temperature, T§ x , the former 
minimum can turn into an absolute one and below this 
temperature, the mean-field mode remains the statistical 
equilibrium, FigJSJd). 



the canonical equilibrium, at T ( 
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FIG. 2. Contour plot of <j> (in log scale), in the (po,pi) space 
for (a) T = 1.4, (b) T = 1.06, (c) T = 0.6 and (d) T = 0.4. 
Simulations realized for Ao = 1 and Ai = 2.3, which yield 



Tf = 1.15, T C = 1 and T C 



0.459. Although p has been 



defined positive, it is also plotted for negative value, for the 
sake of visibility. 



Thus if the apparition of the minimum for the first 
Fourier mode corresponds to a second-order transition 
at T — Tf, then at T = Tq the minimum associated 
with the mean-field mode appears in addition to a saddle 
point between the two minima (spinodal point). Eventu- 
ally, the mean-field minimum turns into an absolute one 
at T = T^, and the canonical equilibrium correspond to 
a new branch. At T = T§ t , the slope of <j) changes during 
the exchange of branch, which is a signature of a first- 
order transition. From a macroscopic point of view, the 
system jumps from a modulated phase to a homogeneous 
magnetized one, as can be seen from the plain lines in 
FigEl In this figure, the local minima of </> are metastable 
states (dotted lines). A mean-field metastable state ap- 
pears when T goes below Tq and the spatially-modulated 
mode turns metastable as the mean-field mode becomes 




FIG. 3. Equilibrium magnetization for the k = (blue) and 
k — 1 (red) modes for the truncated potential: the plain 
lines stand for the equilibrium curves, the dotted ones for the 
metastable states. The symbols correspond to microscopic 
simulations for the complete potential, with N = 1000 and 
Nmc = 10 5 Monte-Carlo steps. 

Molecular simulations confirm these results, as is 
shown by the symbols in Fig|3] Simulations up to 10 8 
Monte-Carlo steps showed that the system either falls 
into the mean-field or the k = 1 phase and is unable to 
escape the potential well it initially falls in. The simula- 
tions were realized with the complete potential, without 
truncation, which confirms that the 2-mode truncation 
approach accepts the main feature of the new configura- 
tion. Modulated phases may also appear in the case of 
long-range magnetization-conserved dynamics 25 1 . 

The metastable states are considered in more detail 
here. The time needed for a long-range system to es- 
cape the metastable state is known to be a function of 
the exponential of N times the barrier in free energy 5<fi 
that the system has to overcome [26} . This barrier here 
corresponds to the difference between the saddle and the 
metastable states. Entropic barriers of order 10 -4 — 10 -3 
were reported in systems of coupled rotators 
here barriers are hundreds times higher, 8(j> 
T = 0.3, making the metastable extremely long-lived. 
Therefore the dynamics of systems with a thousand of 
particles remains trapped in the minimum of free energy 
they initially fall in for very large times, beyond what can 
be achieved today with a desktop computer. FigJ3] illus- 
trates the fact that molecular simuations predict the cor- 
rect amplitude of the equilibrium and metastable states 
(blue and red crosses) up to finite-iV fluctuations, yet 
each realization stayed in the macroscopic state it ini- 
tially falls in, even for long time simulations (10 8 Monte- 
Carlo steps). 

Furthermore, for low temperatures T < 0.5, the two 
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minima are very close compared to the barrier that sepa- 
rate them, so that the typical time for the system to hop 
from the canonical equilibrium to the metastable state 
is of the same order as the time for the reverse process. 
Thus, the canonical equilibrium may be very difficult to 
identify using the time series of order parameters such as 
the population of the mean-field or first Fourier modes. 

The present derivation holds for this Ising-like model, 
yet it can easily be generalized to models with continu- 
ous variables such as systems of rotators with power-law 
potentials 0] 



H = y d-J_Y V{d 



cos 



- Bi), 



(11) 



where (9i,pi) are conjugate variables. Following a similar 
derivation, one can show that the hyperbolic cosine in the 
partition function and in the free energy is replaced by 
the modified Bessel function, Iq [3, H| • Correspondingly 
in Eq.(j8]), the hyperbolic tangent should be replaced by 
Ii/Iq. This set of equations yields the same phenomenol- 
ogy as for the modified Ising model, with the possible 
presence of modulated phases at the canonical equilib- 
rium. The critical temperature of each mode is given by 
T£ = Afe(f + <$o,fc)/4, showing that systems with contin- 
uous variables may as well exhibit modulated phases in 
the presence of competing potentials. 

Finally, considering a more complex case, where 2Ao < 
Ai < ^2fe+i for k' > 1 , three different ordered phases may 
exist at equilibrium and are associated with the mean- 
field phase and to modes f and k' (29[. As before, the 
system switches from one ordered phase to the other via 
a first-order transition. 

Note that the well of free energy the system preferen- 
tially falls in initially, starting from a disordered phase, 
can be deduced from an analysis of linear stability Q 
where each Fourier mode was shown to grow, in the lin- 
ear regime, at an exponential rate that increases with 
Afc. This growth rate, combined to finite- N fluctuations, 
decides in which ordered phase the system initially falls. 
This explains why in Fig|3] more simulations go to the 
modulated phase. Simulations with a lower ratio Ai/Ao 
confirmed that the system distributes more equally be- 
tween each phase during many realizations. 

In conclusion, we have shown that the competition of 
long-range forces leads to new ordered phases, modulated 
at the scale of the system. The canonical equilibrium is 
characterized by first-order phase transition where the 
system switches from one phase to the other, yet the 
different ordered phases also present a strong metastable 
character over such long times that they may be confused 



with the canonical equilibrium. 

The authors acknowledge fruitful discussions with J. 
Barre, G. de Ninno, and especially S. Ruffo. R.B. 
acknowledge support from the Fundacao de Amparo 
Pesquisa do Estado de Sao Paulo (FAPESP). 



[1 
[2 
[3; 
[4 
[5 
[6 
[r 

[8 
[9 

[10 
[11 

[12 

[13 

[14 
[15 
[16 

[17 

[18 
[19 

[20 

[21 

[22 
[23 
[24 
[25 
[26 

[27 
[28 

[29] 



Phys. Rep. 480, 
Phys. Rev. Lett. 80, 5312 
Phys. Rev. Lett. 84, 208 



Campa A., Dauxois T. and Ruffo S. 
57-159 (2009). 

C. Anteneodo and C. Tsallis, 
(1998). 

F. Tamarit and C. Anteneodo, 
(2000). 

A. Campa, A. Giansanti, D. Moroni, Phys. Rev. E 62, 
303 (2000). 

A. Campa, A. Giansanti, and D. Moroni, J. Phys. A 36, 
6897 (2003). 

U. Low, V. J. Emery, K. Fabricius and S. A. Kivelson, 
Phys. Rev. Lett. 72, 1918 (1994). 

S. A. Cannas, P. M. Gleiser and F. A. Tamarit, Recent 
Research Developments in Physics, Vol.5, Part II, pp.751 
(Transworld Research Network, 2004). 
A. Giuliani, J. L. Lebowitz, and E. H. Lieb, Phys. Rev. 
B 74, 064420 (2006). 

A. Giuliani, J. L. Lebowitz, and E. H. Lieb, Phys. Rev. 

B 76, 184426 (2007). 

J.F. Nagle, Phys. Rev. A 2, 2124 (1970). 

D. Mukamel, S. Ruffo and N. Schrieber, Phys. Rev. Lett. 
95, 240604 (2005). 

T. Dauxois, P. de Buyl, L. Lori and S. Ruffo, J. Stat. 
Mech. P06015 (2010). 

E. Edlund and M. Nilsson Jacobi, Phys. Rev. Lett. 105, 
137203 (2010). 

Seul M. and Andelman D., Science 267, 476 (1995). 

G. Malescio and G. Pellicane, Nat. Mater. 2, 97 (2003). 
D. ODell, S. Giovanazzi, G. Kurizki, and V. M. Akulin, 
Phys. Rev. Lett. 84, 5687 (2000). 

A. Dominguez, M. Oettel, and S. Dietrich, Phys. Rev. E 
82, 11402 (2010). 

R. Golestanian, Phys. Rev. Lett. 108, 038303 (2012). 
M. Chalony, J. Barre, B. Marcos, A. Olivetti, and D. 
Wilkowski, larXiv:1202.1258l (2012) . 

Ph.W. Courteille, S. Bux, E. Lucioni, K. Lauber, T. Bi- 
enaim, R. Kaiser, N. Piovella, Euro. Phys. J. D 58, 69 
(2010). 

R. Bachelard, Ph.W. Courteille, R. Kaiser, N. Piovella, 
Europhys. Lett. 97, 14004 (2012). 

F. J. Dyson, Commun. Math. Phys. 12, 91 (1969). 
F. J. Dyson, Commun. Math. Phys. 12, 212 (1969). 
D. J. Thouless, Phys. Rev. 187, 732 (1969) 

T. Mori, Phys. Rev. E 84, 031128 (2011). 

Griffiths R.B., Weng C.Y. and Langer J.S., Phys. Rev. 

149, 301 (1966). 

Antoni M., Ruffo S. and Torcini A., EPL 66, 645 (2004). 
R. Bachelard, T. Dauxois, G. De Ninno, S. Ruffo and F. 
Staniscia, Phys. Rev. E 83, 061132 (2011). 
Note that since the tanh is an odd function, Eq. ((8} leaves 
the odd and even modes uncoupled. 



